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Abstract. Various bounds for the energy of collective excitations in the Heisenberg 
antiferromagnet are presented and discussed using the formalism of sum rules. We 
show that the Feynman approximation significantly overestimates (by about 30% in 
the S = ^ square lattice) the spin velocity due to the non negligible contribution 
of multi magnons to the energy weighted sum rule. We also discuss a different, 
Goldstone type bound depending explicitly on the order parameter (staggered mag- 
netization). This bound is shown to be proportional to the dispersion of classical 
spin wave theory with a q-independent normalization factor. Rigorous bounds for 
the excitation energies in the anisotropic Heisenberg model are also presented. 
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1. INTRODUCTION 



In the last few years a considerable number of papers has been devoted to the 
study of the Heisenberg model for antiferromagnetism, especially in 2-D. This inter- 
est is mainly motivated by the need for a better understanding of the antiferromag- 
netic behavior of the undoped precursor insulators of the high superconductors. 
After the pioneering works by Anderson [1] and Oguchi [2] , based on spin wave the- 
ory, several theoretical methods have been developed to study this problem. These 
range from spin-wave theory up to second order in 2^ to series expansion methods 
from the Ising side and to Monte Carlo calculations (see the review papers [3-4] for 
exaustive discussions and references). 

The purpose of this paper is to discuss the elementary excitations of the Heisen- 
berg antiferromagnet at zero temperature using a sum rule approach. Only recently 
systematic theoretical investigations of the dispersion of spin waves in the whole 
Brillouin zone have become available [5-10]. Recent experiments [11] in La2Cu04 
with neutron scattering suggest that the dispersion follows the predictions of classi- 
cal spin wave theory with a proper renormalization factor. Even at low q, where the 
dispersion becomes linear, rather relevant questions still remain to be clarified in a 
satisfying way. Among them we recall the problem of the validity of the so called 
"Feynman" or single mode approximation for the calculation of the spin velocity 
and of the role of multi magnon excitations. These questions, first discussed by 
Hohenberg and Brinkman many years ago in ID antiferromagnets [12], have been 
recently addressed by Singh [9] in the S = 1/2 square lattice. In the present work 
we are mainly interested in the 2D case and in general in systems with broken 
symmetries. 

The paper is organized as follows: in sect. 2 we discuss the Feynman approach 
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to spin excitations and we prove that it cannot reproduce the correct dispersion 
of spin waves at low q because of the presence of multi-magnon excitations which 
affect the energy weighted sum rule also in the long wave length limit. In sect. 3 we 
discuss a different bound for the energy of elementary excitations. This bound, first 
introduced by Wagner [13] many years ago, has the form of a Goldstone theorem and 
depends explicitly on the order parameter. It can be easily calculated through the 
whole Brillouin zone and in particular it exhibits the same dependence on q as the 
one given by classical spin wave theory (SWT), with a proper renormalization factor. 
In sect. 4 we present results for the anysotropic Heisenberg model. In particular 
we derive rigorous upper bounds for the mass gap in the easy-axis antiferromagnet 
and for the gapless dispersion law in the case of the easy- plane antiferromagnet. 

2. THE FEYNMAN APPROXIMATION 

In the following we investigate spin excitations in the framework of the Heisen- 
berg model for anti-ferromagnetism (AFM) characterized by the Hamiltonian 

<ij> 

where < ij > denotes a sum over all nearest-neighbor pairs and J > 0. The limits 
A = and A = 1 correspond to the most famous Ising and isotropic Heisenberg 
models respectively. At zero temperature the isotropic Heisenberg model is believed 
to give rise to spontaneous sublattice magnetizazion also in 2D (square lattice), 
though quantum fluctuations have a crucial role in reducing the value of the order 
parameter (actually the S = ^ Heisenberg model has been rigorously proven to give 
rise to spontaneous magnetization only in 3D [14]). In sections 2. and 3. we mainly 
discuss the isotropic case (A = 1) and we assume the staggered magnetization to 
be oriented along the z-axis. This is also the case for the anisotropic case if A < 1. 
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Conversely when A > 1 (see sect. 4.) the axis of (spontaneous) magnetization hes 
in the x — y plane (easy plane). 

In the following we will mainly consider excitations generated by the spin op- 
erator: 

< = ;^E<^^"'" (2) 

^ i 

These excitations are transverse with respect to the ^-staggered magnetization axis. 
The most important among such excitations are spin waves (magnons) that repre- 
sent the elementary excitations of the system. Rigorous upper bounds for the energy 
of these excitations can be obtained at zero temperature using the sum rule method. 

The most popular bound is given by the Bijl-Feynman ansatz, analog of the 
most famous approach employed to investigate the propagation of density excita- 
tions in Bose superfluids [15]. It is obtained by applying the spin operator (2) to 
the ground state of the system: One finds: 

I F >= > ^ I > (3) 

In eq.(3) S^{q) =< | s^q, Sq | > is the transverse structure factor entering here 
as a normalization factor. The excitation energy of the "Feynman" state is given 

by 

1< I [«-a. [H:S^]] I > 

e^iq) =<F\H\F>-<0\H\0>= g^q) 

and provides, at zero temperature, a rigorous upper bound for the energy e{q) 
of the lowest state excited by the operator Sq. This can be directly shown by 
identifying the numerator and the denominator of eq. (4) as the energy- weighted and 
non energy- weighted moments of the transverse dynamic structure factor {q,i^) = 
|<0|Sq|n>p5(a; — a;„o)- In fact, using the completeness relation, one can 
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write 



J S^{q,u)udw = ^ |< I I n >|2 cu^o = ^ < I [s^^, | > (5) 



and 



j S^{q,ij)diJ = ^ |< I s^q I n >|2=< I s^qS^ I >= S^{q) (6) 



Note that at T = the dynamic structure factor S-^{q,Lij) vanishes for a; < 0. 

The Feynman energy (4) has been aheady used by several authors to study the 
energy of elementary excitations in the Heisenberg model [12,5-6,9]. The numerator 
of eq.(4) can be easily calculated employing the commutation rules for the spin 
operators. The result is 
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< I [s\, [H, s^]] I >= z[Ml - A7q) + /,(A - 7q)] (7a) 



Analogously, for the and operators one finds: 

^ < I [sl^, [H, sy^]] I >= z[Ml - A7q) + m - 7q)] (76) 

^ < I [slq, [H, 4]] I >= Xzif, + fy){l - 7q) (7c) 

where z is the number of nearest neighbors, 7q = ^ Yls cosq • 6 and we have intro- 
duced the quantities 



fx — 


*f*H-5 > 


/.= 




fz = - 


~"2 < > 



(8) 



Here S is the lattice vector connecting nearest neighbors. In the square lattice one 
has 7q = ^{cosqx + cosqy) while in the cubic lattice 7q = \^{cosqx + cosqy -\- cosq^)-, 
having set the lattice parameter equal to 1. It is worth noticing that the form of 



the energy weighted sum rule relative to Sq differs from the one relative to Sq and 
s^. This follows from the fact that the Heisenberg Hamiltonian (1) is invariant for 
spin rotation in the x — y plane. 

In the isotropic case (A = 1) eq.(7a) becomes 

^ < I [.^q, [H, I >= Z{f, + fy){l - 7q) . (9) 

Note that even in the isotropic limit A = 1 the quantity fg differs from fy{= fx) if 
there is spontaneous magnetization along the ^-axis. 

At small q the energy weighted sum rule (9) becomes (we consider here for 
simplicity the square and cubic lattices where 7q = 1 + ^q-^ + 0(g'^)) ): 

^<0|[.^q,[i?,<]]|0>=(/, + /,)5^ (10) 

and exhibits the typical dependence characterizing the most famous /-sum rule 
for density excitations [16]. 

The denominator of eq. (4) is the Fourier transform of the two-body transverse 
spin correlation function. Its behavior is dominated, at low q, by long range correla- 
tions associated with spin waves. Numerical results for S-^{q), based on Monte Carlo 
calculations [5-6] and series expansion methods [9], are now becoming available. 

From a general point of view the Feynman energy (4) is expected to provide 
a good estimate for the frequency of elementary excitations in Heisenberg antifer- 
romagnets. This system can be in fact considered a relatively weakly interacting 
many body system as compared, for example, to other strongly interacting quan- 
tum system such as superfluid '^He where the Feynman approximation is known 
to overestimate in a significant way the the energy of lowest excitations at high 
momenta. 

An important question is howevever to understand what happens to the Feyn- 
man approximation in the long wave length limit dominated by the propagation 
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of macroscopic spin waves. While in superfluid ^He the Feynman ansatz is known 
to reproduce exactly the phonon dispersion (in terms of sum rules this means that 
both the energy weighted and non energy weighted sum rules for the density oper- 
ator are exhausted by phonons) the situation is different in the spin case. In fact 
the non conservation of the spin current makes the contribution of multi-magnon 
excitations particularly important in the low q limit. These excitations exhaust 
a finite fraction of the energy weighted sum rule (EWSR) and consequently the 
Feynman energy (4) does not approach the correct dispersion law at small q. In 
the following we will discuss such an effect in a quantitative way with the help of 
available microscopic calculations of the spin stiffness coefficient. 

It is convenient to write the transverse dynamic spin structure function in the 
following way 

(g, CO) = Aiq)diu; - iviq)) + S^^iq, co) (11) 

where we have separated the sharply peaked single magnon contribution character- 
ized by the dipersion law u>{q) and strength A{q), from the smooth contribution 
S;^^{q,u}) arising from multi magnon excitations (S'^^(g, a;) = for a; < a;(g)). 

The main results for the single magnon and multi magnon contributions to the 
various moments of S-^{q, co) at small q are summarized in table 1. The main point is 
the q^ dependence of the strength associated with multi-magnon excitations. This 
dependence differs from the q'^ dependence associated, for example, with multi- 
phonon excitations in Bose superfluids. The difference is due to the fact that the 
current is conseved in Bose superfluids because of translational invariance. In the 
case of spin excitations the quantity [H, s^] , proportional to the spin current (see 
eq.(15) below), is not conserved even in the low q-limit and this implies a stronger 
q'^ dependence for the strength associated with multi-magnon excitations. A similar 
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behavior is exhibited by spin excitations in normal Fermi hquids [17]. This result 
implies that multi-magnon excitations affect the energy-weighted sum rule with a 
term proportional to [18]. 

The occurrence of a contribution to the energy weighted sum rule due to 
multi-magnon excitations is clearly exploited by the calculation of the double com- 
mutator relative to the "longitudinal" operator Sq = Si -^f s*'* (see eq.(7c)) 
for which we find, at low q, 

^<o|[.iq,[iy,<]]|o>=(/, + /,)5^ (12) 

This contribution, quadratic in q, is entirely fixed by multi-magnon excitations since 
single magnons are not excited by Sq. 

The low q contribution to the tranverse energy weighted sum rule (5,10) arising 
from single magnons is given by \psq^ where ps is the spin stiffness coefficient. This 
can be easily understood by using the hydrodynamic expression for the spin velocity 
[19]: 

where 

Xdq) = 2j2\<0\s^_jn>\'^ = 2 [du.^^^ (14) 

is the transverse magnetic susceptibility. This sum rule is expected to be entirely 
exhausted, at low q, by the one magnon excitation. If the energy weighted sum 
rule (5,10) were also entirely exhausted by the one magnon mode at low q, then the 
ratio 

lirrin 



1 En 


l< 


1 «-q 


n > 


|2 , , 




< 1 


S-d 1 


n > p /ijJnO 



X^(0) 

should coincide with c^. The comparison between the quantities 2(/z + fy) and 
Ps then provides a direct and quantitative information about the contribution of 
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multi magnons to the energy weighted sum rule. Both the quantities (/z + fy) 
and ps are now available through different theoretical calculations. All the various 
predictions, based on spin wave theory to second order in 2^ [7,20], series expansion 
from the Ising side [21] and Monte Carlo calculations [5-6] agree with the value 
2{fz + fy) = 0.25 in the S = ^ square lattice. Viceversa the most recent estimates 
for Ps [20-22] predict values in the range 0.18 — 0.20. Since the non energy weighted 
sum rule (6), entering the denominator of the Feynman bound (4), is expected 
to be exhausted by the single magnon (see eq.(19) below), we then conclude that 
the Feynman ansatz overestimates the spin velocity by about about 30%. In the 
S = 1 square lattice the overestimate is about 10%. In table 2 we report, for 
completeness, the values of various thermodynamic parameters relative to the 2D 
Heisenberg model. These values correspond to the predictions of spin wave theory 
up to ^2g)2 [23] and are rather close to the ones given by the series expansion method 
from the Ising side and by Monte Carlo calculations. 

It is useful to study more explicitly the role of the spin current and its connec- 
tion with the spin stiffness coefficient and the energy weighted sum rule. To this 
aim let us start from the continuity equation for the spin density (in the folUowing 
the vector q will be taken along the x-axis): 

<ij> 

definining the component of the spin current parallel to q. Equation (15) provides 
the following expression for the spin current at q = 0. 

if.(o) = -;4E<^r+A (16) 

i,S 

where = Xi — Xj is the x component of the vector connecting the nearest-neighbor 
pair < ji >. 
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The key point is that the spin current (16) is not a conserved quantity (it does 
nor commute with the Hamiltonian) and consequently, when apphed to the ground 
state, it can give rise to excitations with non-vanishing strength. Such excitations 
are multi-magnon states since spin waves with q = cannot propagate. 

Let us now calculate the static response relative to the current jg^{q). Due to 
the equation of continuity (15), this is exactly fixed by the energy weighted sum 
rule for the spin operator 

= 4 5^ l< |< I n >|2 a;„o = 2(/, + fy) (17) 

^ n 

where we have taken the low q limit (10) of the energy weighted sum rule. Both spin 
waves and multi-magnon excitations affect this quantity at low q. The spin wave 
contribution is fixed by the spin stiffness coefficient (see the discussion above and 
table 1), while the multi-magnon contribution can be calculated through the static 
response of the q = component (16) of the spin current operator. In conclusion 
we get 

Ps = 2{h + fy)-x{jim (18) 

Result (18) for the spin stiffness coefficient ps shares important analogies with the 
most famous expression Ps = P — Pn for the superfluid density of a Bose liquid. In 
eq.(18) the quantity 2{fz + fy) plays the role of the total density p, fixed by the 
model independent /-sum rule [16], while the quantity xUs^i^)) plays the role of 
the normal density pn, defined as the low q limit of the transverse current reponse 
function [24]. Note that in the case of antiferromagnetism, where the current is 
not conserved, we can safely take the q ^ limit of the current operator for the 
calculation of the multi magnon contribution to the static current response. 
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It is remarkable to point out that relation (18) was obtained in an independent 
way by Singh and Huse [21] starting directly from the definition of the spin stiff- 
ness as helicity modulus. The full agreement between the two formal derivations 
provides further support to the theory of spin hydrodynamics and at the same time 
emphasizes the role played by multi magnon excitations. Concerning this last point 
it is worth noting that in the large S limit multi magnon excitations are absent, 
X(jf^(0)) = and ps coincides with 2{fg + fy). Actually, using the results of spin 
wave theory [23], one can easily show that the multi magnon term xUs^i^)) is sec- 
ond order in while the longitudinal sum rule (12), dominated by multi magnons, 
is first order in This different behavior is likely associated with the fact that 
longitudinal excitations are mainly two magnon states, while the multi magnon 
component of the transverse response is dominated by three magnon states. 

Another important result emerging from table 1 concerns the low q behavior 
of the transverse spin structure factor (6): 

S^iqh^O = l^q (19) 

accounting for the fluctuations associated with the propagation of long wavelength 
spin waves. The coefficient of linearity has been directly calculated by Singh [9] using 
the series expansion method. The resulting estimate is in reasonable agreement with 
eq.(19). 

It is finally useful to stress that the results discussed in this section using the 
sum rule technique emphasize in an explicit way the existence of a spontaneously 
broken symmetry in spin space. Different results would be obtained if one instead 
decided to work with an isotropic ground state, as happens, for example, in a 
numerical simulation in a finite system. In this case the results for the excitation 
energies, obtained through the evaluation of sum rules, would correspond to an 
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average between transverse and longitudinal excitations and the information on the 
dispersion law of elementary modes would be consequently poorer. 

3. ORDER PARAMETER AND EXCITATION ENERGIES 

The discussion of sect. 2. on the behavior of the Feynman energy in the low 
q region is based on the analysis of the spin structure function. The existence 
of spin waves with linear dispersion must be however assumed in order to discuss 
such a behavior and cannot be predicted using this method, unless one exploits 
numerically the rather difficult low q-regime. For this reason it is useful to derive 
alternative bounds for the excitation energies which exploit more directly the low 
q regime. Such bounds can be obtained with the help of an inequality due to 
Bogoliubov and point out a crucial feature characterizing antiferromagnets as well 
as other systems with spontaneously broken symmetries: the existence of an order 
parameter. This phenomenon is known to be at the origin of Goldstone modes 
which, in the antiferromagnetic case, take the form of spin waves with a linear 
dispersion at low q. This approach was first proposed by Wagner [13] to prove 
the existence of Goldstone modes in an important class of physical systems. To 
our knowledge it has never been used to investigate the full q-dependence of the 
excitation spectrum of Heisenberg antiferromagnets. 

The starting point is the introduction of an upper bound for the energy u;{q) of 
the lowest excitation with wave vector q, in terms of the ratio between the energy 
weighted and the inverse energy weighted sum rules relative to the operator s^: 



In eq.(20) we have made use of eq.(5) and used definition (14) for the transverse 
susceptibility. 




(20) 
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The upper bound (20), holding at zero temperature, is stronger than the Feyn- 
man one (see eq.(4)), being based on the inverse energy weighted sum rule X'^il) 
rather than on the non energy weighted sum rule S^{q). Its determination requires 
however the difficult calculation of the q-dependence %"'"(?)• the following we 
will combine the bound (20) with the Bogoliubov inequality [13,25] for the static 
response relative to the operator 

X^iq) < >>l< [^-q.<-g] >\" (21) 

This inequality introduces the "conjugate" operator Sq_g where g is the antiferro- 
magnetic vector fixed by the condition e*^'^ = 1 when R connects sites in the same 
sublattice and —1 when it connects sites in different sublattices. 

Using inequality (20) and (21) we then obtain the useful rigorous result [26] 

A major advantage of inequality (22) as compared to the Feynman bound (4), is 
that it involves commutators both in the numerator and denominator. In particular 
the quantity 

< [^-q, >= ^ < ^ E ^"^ (23) 

i 

coincides with the staggered magnetization (assumed here along the 2;-axis), i.e. 
with the order parameter of the problem, and is independent of q. 

The full q'-dependence of the bound (22) is then entirely fixed by the double 
commutators entering the numerator. Such commutators have been already calcu- 
lated in sect. 2 (see eq.(7)). Noting that 7q-g = — 7q we find the following result 

< ^/ir^ = !(4±A)<,s„.(,) (24) 
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where (jj^^{q) = zJS^Jl — 7q is the dispersion law of classical spin wave (SW) 
theory [1] and we have used the property fx = fy 
The following remarks are in order here: 

i) The rigorous bound (24) exhibits a linear behavior in for — > 0, provided 
the order parameter is different from zero (Goldstone theorem). Furthermore this 
bound is symmetric by exchange of q with g — q and hence predicts the vanishing 
of elementary excitations also at the staggered wave vector g. 

ii) The q-dependence of this bound is entirely contained in the classical law 
uj^^{q), the coefficient of proportionality being independent of q. In particular 
from eq.(24) we obtain the bound 

c < ^iklMcSW = 2V2^^^^^^tA) (25) 
SmJ m 

for the spin velocity in terms of the quantities {fz + fy) and m {c^^ = y/2zSJ is 
the prediction of classical SW theory). Using the numerical results of table 1 for 
{fz + fy) and m we find c < 1.60"^^ in the S — ^ square lattice. The bound (25) 
overestimates by ~ 30 % the value of the spin velovity calculated through equation 
(13) (c = 1.2c^^). In the S = 1 square lattice result (25) yields c < 1.2c^^ while 
eq.(13) gives c = l.lc"^^. At small q the quality of the new bound is hence similar 
to the one of the Feynman approximation. From a conceptual point of view it has 
the advantage of exploiting directly the low q behavior with the only assumption of 
the existence of a broken symmetry. It is also interesting to remark that, using the 
result of second order spin wave theory [23], the bound (25) for the spin velocity 
coincides with the exact value (13) up to first order terms in Deviations from 
the exact value are associated with multi magnon effects (terms in ^^^^ ). 

The dispersion of magnon excitations in the S = ^ square lattice Heisenberg 
model has been the object of a recent Monte Carlo calculation [5]. The authors 
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of ref.[8] have fitted their results with the law a;(g) ~ 1.2(jj^^{q) (similar results 
have been very recently found also by the authors of ref.[10]), consistently with the 
value of the spin velocity obtained from eq.(13). The upper bound (24) is then 
found to overestimate the magnon dispersion by the same amount (~ 30 %) in 
the whole Brillouin zone. In fig.l we report the prediction of the Goldstone-type 
bound (24) together with the fit to the results of ref. [8] and the predictions of the 
the Feynman approximation taken from ref. [9]. It is interesting to remark that the 
Feynman approximation is much more accurate near the maximum of the dispersion 
curve rather than in the low q region where, according to the discussion of sect. 2, 
it overestimates the linear dispersion by ~ 30 %. 

iii) Inequality (24) becomes an identity in the large S limit (/^ = ^S'^, fx = 
fy — 0, m — S) where it coincides with the prediction classical spin wave theory 

[!]• 

The Bogoliubov inequality (21) can be used to provide directly a bound for the 
transverse susceptibility x^{q)- Using the relation 7q-g = — 7q one finds 



2 

m 



At 5 = eq.(26) yields 



^'^"^^-MITTM' '''' 



2 



while near the staggered vector g one finds the typical divergent behavior 

characterizing the transverse staggered susceptibility. 

Once more these inequality become identities if one works with spin wave theory 
up to first order in ^ . Deviations from the exact results for these formulae are the 
direct consequence of the role of multi magnon excitations. 
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It is finally useful to complete the analysis of sect. 2 concerning the contribution 
to the various sum rules given by the single magnon and multi magnon excitations 
in the region of the staggered vector g. The results are reported in table 1. We 
note that single magnons exhaust the transverse structure factor and susceptibility 
sum rules characterized by typical infrared divergencies. The result for the spin 
structure factor near the staggered vector can be obtained with the help of the sum 
rule (23) 

y^[< I s^q \ n><n \ Sq_g I > - < I s\_^ | n >< n | s!iq | >] 

n 

=<[s%,sl_^]>=im (29) 

In fact, since the magnon matrix element < | \ n > behaves like y/q at low 
q (see table 1 and eq.(19)), it follows that the sum rule (29) can be satisfied only by 
a divergent behavior of the magnon matrix element < n \ Sq_g | > (multi magnon 
excitations give rise to higher order contributions) according to the equation 



< n I I >=< I I n' >= \ ^ (30) 

2<0| s^q I n > 

holding for g — > 0. Here | n > and | n' > are single magnon states with opposite 
wave vector and we have assumed, without any loss of generality, the matrix element 
< I s^liq I n >=< n' I I > to be real. The magnon contribution (30) 
dominates the divergent behavior of the spin structure factor near the staggered 
vector that then takes the form: 



5^(|g-q|).^o = |^ (31) 



The above results are consistent with the rigorous inequality [27] 

S^{q)S^{\g-ci\)>lm' (32) 
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following from the uncertainty principle and holding for any value of q and for 
any antiferromagnetic system. According to results (19) and (31), the uncertainty 
principle inequality becomes an identity in the q ^ limit. The coefficient of the 
^ law (31) has been recently calculated in the S = ^ square lattice by Singh [9] 
using using the series expansion method from the Ising side. His prediction turns 
out to be larger by (~ 20%) than the value predicted by eq.(31). This discrepancy 
remains to be understood. 

Result (31) can be used to study the quality of the Feynman energy (4) near 
the staggered vector g. One finds: 

'^Fd g-q Dq^o = 1 cq (33) 

where we have used expression (13) for the spin velocity c. Result (33) overestimate 
the spin velocity by ~ 30 % in the -S' = ^ square lattice. The enhancent coincides 
with the ratio between the left and right hand sides of inequality (27) for the trans- 
verse suscptibility and follows from the multi magnon contribution to the energy 
weighted sum rule. 
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4. RESULTS FOR THE ANISOTROPIC HEISENBERG MEDEL 



The energy weighted sum rule (7) for the Heisenberg model has an interesting 
behavior at low q in the anisotropic case (A 7^ 1). In fact at q = eqs.(7a) and (7b) 
become: 

hm i < I [s^q, [H, s^]] I >= z{l - A)(/, - /,) (34a) 



and 



hm l<0\ [sl^, [H, s^]] I >= z{l - A)(/, - U) (346) 

Conversely the EWSR relative to vanishes with q since the Heisenberg Hamilto- 
nian (1) conserves the 2;-component of the spin operator. 

Note that the quantities fz — fy and fz — fx must be positive for A < 1 and 
negative for A > 1. This is a rigorous stability criterium imposed by the positivity 
of the energy weighted sum rules (34). 

Result (34) can be used to derive a rigorous upper bound for the mass gap 
when A < 1. In fact in this case eq.(22) yields 



2z 



a;(g = 0)<— ./(/|-/2)Vl-A2 (35) 



where we have explicitly specified that the magnetization is along the 2;-axis (easy 
axis) and used the property 7g —1. 



This upper bound exhibits the typical non analytic Vl — A^ behavior predicted 
by SWT near A = 1. In the S = ^ square lattice the coefficient of proportionality of 
the upper bound (35) is equal to 1.9, compared to the value 1.3 obtained in ref.[28] 
using the series expansion method. 

Using the Bogoliubov inequality (21) it is also possible to obtain the rigorous 
bound 
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for the transverse staggered suceptibility. 

Both results (35) and (36) apply only to the case A < 1. It is also interesting 
to discuss the behavior of the system beyond the isotropic point A = 1 where one 
expects the spontaneous magnetization to occur in the x — y plane (easy plane). In 
the following we assume the magnetization axis to coincide with the a;-axis. One 
can find in this case a rigorous Goldstone type upper bound similar to eq.(24). This 
bound is obtained starting from inequality (22), by replacing the operator Sg with 
Sq (the replacement follows from the new direction of the magnetization axis): 

... < <i^i,,ig,<]]><K-,.i-ff,<-j> 
" i<[«i„<-j>p '''^ 

Using results (7) for the corresponding double commutators and the identity 
< ["Siq, ■Sq-g] ~^ < ;^ Si -sf e*^ '"* >= —inix (staggered magnetization along 
the X-axis), we obtain 

uj\q) < ^(/. + /,)(! - 7q)[/.(l + A7q) + /.(A + 7q)] (38) . 

yielding a linear dispersion for uj{q) at small q (the occurrence of gapless spin exci- 
tations for the easy plane antiferromagnet has been recently pointed out in ref.[29]). 
It is worth noticing however that, differently from eq.(24) holding in the isotropic 
case, the bound (38) is not symmetric by change of q with g — q and in particular 
it is not gapless at the staggered point g. This reflects the fact that this system, 
characterized by an anisotropy of the Hamiltonian in the ^-direction and by a spon- 
taneous staggered magnetization along the a;-axis, exhibits two different branches in 
the excitation spectrum: one excited by the operator Sq and for which eqs. (37-38) 
provide a rigorous upper bound, and one excited by the operator s^. The bound 
for the second branch is easily obtained by replacing, in eq.(37), the operator Sq 
with and Sg_q with Sg_q. This corresponds to replacing q with g — q and hence, 
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in eq.(38), 7q with — 7q. Notice that this second branch is gapless at the staggered 
vector g. 

Equation (38) provides a rigorous upper bound for the spin velocity holding 
for an arbitrary value of A (larger than 1 of course): 



c < V^A(1 + A)J(/, + /,)(/, + /,) (39) 



Result (39) coincides with result (26) in the A — > 1 limit and provides a non trivial 
result also in the A ^ oo limit (XY model). 

Another interesting result can be obtained for the behavior of the derivative 
of the energy with respect to the transverse coupling constant A. This behavior is 
important because it characterizes the nature of the phase transition. The derivative 
can be calculated starting from the general Feynman formula 

dE{X) 



dX 



-Z{f. + fy) (40) 



which straightforwardly follows from the form of the Heisenberg Hamiltonian (1) 
and definitions (8) for fx and fy. When A ^ 1~ one has fx=fy¥'f7i while 
when A — > 1"*" one has /+ = /~ and ft — fy — fy ■ This finally yields 

dE{\y 



dX 

dE{X) + 



-'2zf, 



= -<f7+fy)- (41) 



dX 

Using the values for and fy reported in table 2 (corresponding to spontaneous 
magnetization along the 2;-axis and hence to /~ and /~ respectively) we find 
^^^^ = —0.32 and '^^^^^ = —0.50. These values are in excellent agreement 
with the results obtained in ref.[30] through a direct Monte Carlo calculation of the 
energy as a function of the coupling constant A. 
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CONCLUSIONS 



In the present work we have derived several new results concerning the propa- 
gation of elementary excitations in the Heisenberg antiferromagnet. In particular: 

1) We have proven that the Feynman approximation does not yield the correct 
dispersion of long wavelength spin waves, due to the role of multi magnon excita- 
tions which contribute to the energy weighted sum rule (EWSR) even in the low q 
limit. Physically this behavior originates from the fact that the spin current is not 
conserved. Actually the multi-magnon contribution to the EWSR is fixed by the 
static spin current polarizability x(jf^(0)) (see eq.(18)). Due to this effect, second 
order in 2^ , the Feynman approximation turns out to overestimate the spin velocity 
in the S = \ square lattice by about 30%. 

2) We have derived (sect. 3) a Goldstone-type bound for the energy of spin exci- 
tations. This rigorous bound depends explicitly on the order parameter (staggered 
magnetization) and is proportional to the classical dispersion of spin wave theory 
with a q independent normalization factor. It consequently vanishes at q = as well 
as at the staggered wave vector q = g. This bound is shown to have an accuracy 
similar to the one of the Feynman approximation. 

3) We have obtained useful results also for the anisotropic case (sect. 4). In 
particular for the easy-axis antiferromagnet we have derived a rigorous bound for 
the mass gap. Viceversa the upper bound in the easy plane antiferromagnet is 
proven to be gapless in agreement with the general statement of the Goldstone 
theorem. We have also explicitly calculated the discontinuity of the derivative of 
the energy with respect to the transverse coupling constant at the isotropic point. 

A more systematic investigation of the structure of elementary excitations in 
the anisotropic case (including the X-Y model) will be presented in a future paper. 
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FIGURE CAPTION 

Dispersion of spin excitations in the S = ^ square lattice {q^ = %)■ The long- 
dashed line corresponds to the ijj{q) = 1.2(jj^^{q) fit to the Monte Carlo results of 
ref.[8,10]; the squares (taken from ref.[9]) correspond the Feynman bound (4), while 
the full line to the Goldstone-type bound (24). The prediction of classical spin wave 
theory Lij{q) = Lij^^{q) is also repoted (dashed line). 
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TABLE 1 







magnon 


multi-magnons 








const 


IK 


\ |2 
)nO 1 






1(4- 


-qJnO 1 


2cm?' / PsQ 


const 


T 1 


(■Sq)nO P /"^nO 




q' 


T 1 


(Sq)nO P 


paq/2c 


q' 




(^q)nO P <^no 




q' 


T 1 


(Sg_q)nO P /^nO 




const 


Sn 1 


(Sg_q)nO P 


m^c/ipsq 


const 




(Sg_q)nO P <^nO 


m^c^ jlpa 


const 



Matrix elements, excitation energies and sum rule contributions from one-magnon 
and multi-magnon excitations at T = 0. 
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TABLE 2 







E 


m 




Ps 


c 


fz + fy 


fz fy 


s = 


1 

2 


-0.67 


0.30 


0.061 


0.18 


1.7 


0.125 


0.04 


s = 


1 


-2.33 


0.80 


0.092 


0.87 


3.1 


0.47 


0.25 



Parameters of the isotropic 2D AF Heisenberg model predicted by spin wave theory 
up to second order in (25)^ [^^]- '^^^ Heisenberg couphng constant J has been set 
equal to 1 and magnetization is taken along the z-axis. 
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